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Abstract: 

A large CP-asymmetry AAqp has been reported in the D° — > 7t + tt~ /K + K~ system. At 
present it remains unclear whether this is due to incalculable strong interaction matrix 
elements or genuine new physics (NP). Amongst the latter a new weak phase in the 
chromomagnetic operator Os has emerged as a promising candidate. Extending earlier 
ideas we show that the interference of long-distance (LD) terms with the Og matrix 
element, which has a large strong phase, gives rise to direct CP-violation at the level of 
a few percent in D° — > (p°,u) , ~f and Dj, — > (p + ,K* + )'~f for reference values lm[C$ p ] ~ 
0.4 • 10 -2 . This is two orders of magnitude above a Standard Model (SM) estimate. The 
contribution of Im[C^ p ], which is dependent on the model of NP, is governed by the 
LD strong phase which vanishes in the chiral limit at leading order. The question of 
whether this is significantly changed by radiative corrections is an open and interesting 
question that we discuss. Furthermore we point out that the relative size of left- and 
right-handed (photon polarisation)-LD amplitudes can be measured, in principle, through 
time- dependent CP (TDCP) asymmetries in the case where they are both sizeable which 
is supported by SM estimates. Thus determination of the latter provides interesting 
information on the LD-chirality independent of NP. We comment on the origin of the 
LD contribution, which we believe to be dominated by weak annihilation (WA), in the 
appendix. 
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1 Introduction 



1.1 A cv in D° tttt/KK 

Recent LHCb [2] and CDF [3] data both report significant evidence for CP-violation in 
the charm system. The combined value leads to a world average [3] 

AA CP = AS K ~ - a£?~ = -0.65(18) ■ 1(T 2 , (1) 

where 

A t __ nD o ^f]-np°^f] (2) 

is a shorthand for the time integrated CP-asymmetry, for a case where the final state / 
is a CP-eigenstate. AAcp is a convenient quantity since systematic experimental errors 
cancel. It is worthwhile to add that if SU(3)p, or more precisely [/-spin, were a good 
symmetry then A^p K = — A 1 ^ . In the quantity the AAqp TDCP-asymmetry part 
cancels. Effects can though remain through time-acceptance differences in the tt- and 
.fT-system. The latter is estimated to be small, e.g. [2], and thus in summary direct ( i.e. 
time-independent CP-asymmetry) is expected to be responsible for the relatively large 
value of AAqp- 

Sizeable direct CP-asymmetries, c.f. appendix [Cj necessitate large strong (CP-even) 
and weak (CP-odd) phase differences in two amplitudes of comparable size. The rea- 
son CP-violation is believed to small in the charm system is that the weak phases 
are suppressed by four powers of the Cabibbo angle, leading to the naive expectation 
AAqp ~ few • 10 -4 . In the non-leptonic case the QCD matrix elements, which determine 
the strong phase as well as the ratio of amplitudes, are difficult to compute from first 
principles as the size of the charm mass is neither suited to chiral nor heavy quark the- 
ory. Thus the question of whether the large central value (JI]), should it remain, is due 
to NP jU [5j El [7] or somewhat unexpected strong dynamics [8j |9j [TOj [11] , such as in the 
AI = 1/2-rule K — > txtx systerrj^j is an open question at present. It is fortunate that 
contributions due to new AI = 3/2-operators can be tested against the SM with isospin 
sum rules [12] . 

Taking the viewpoint that the asymmetry is largely due to NP it turns out that a 

x It was pointed out quite some time ago [14] that an enhancement of the triplet transition, in the 
SU(3)-flavour classification, may lead to sizeable CP-violation. E.g. A^p ~ 0.08 • 10~ 2 which would lead 
to |AAcp| — 0.16 • 10~ 2 which is not far off the central value in ([T]). 
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weak phase in the |AC| = 1 chromomagnetic operatoij^] 



8 = -^ua.G(l+ l5 )c, 0' 8 = -^ua.G(l- l5 )c (3) 



(a ■ G — a flu G^X a /2), appears to be a promising candidate [15]. By which we mean 
that it does not contradict observations elsewhere, such as D°-D° -mixing. Note, the O s ^- 
operators are of the AI = 1/2-type and do not fall into the testable AI = 3/2-class 
mentioned earlier. Furthermore, 0' s is the structure which is the less abundant helicity in 
the SM due left-handedness of the weak interactions; [Cg/CsjlsM ~ m u /m c . 

To get an idea of the size of the NP contribution [16J one might resort to naive 
factorisation (NF) e.g. [13J. Slightly extending the notation in [TB] one gets, 

AA™p | NF ~ -1.8 (lm[Cf p ] - lm[Cf s NP ]j sin(S) , (4) 

where 5 is the unknown strong phase difference between the KK and tttt states which is 
expected to be sizeable. Note since the sign of sin S is unknown the sign of the difference 
of the two Wilson coefficients is currently not determined by the D° — > ttti/ KK-system. 
Throughout this paper, if not otherwise stated, the Wilson coefficient are understood to 
be evaluated at the charm scale. Since the decay of a J P (D°) = CT particle into two 
J p (n/K) = CT particles necessitates parity violation only the 75-part in ^ contributes 
and therefore results in opposite signs of Im[Cg^ p ] and ImfCg^] in Q respectively. It is 
noted that a value of 

flm[Cf p ] - lm[C 8 NP ]\ sin(5) ~ 0.3 ■ 1(T 2 naive factorisation (NF) (5) 

could account for the central number in Q. One has to bear in mind that ^ is due 
to NF and could easily be out by a few factors. Following the literature we shall take 
lm[CP NP ] = 0.4 • 10- 2 reference value, which is two to three orders of magnitude 



above the SM-value for C 8 , c.f. appendix |B.1[ and suppressed by an additional factor 
m u /m c in the case of C 8 . 



1.2 ^ C p in D° I/7 

As lamented above the situation in D° — > mi / KK remains unclear. Thus the question of 
whether a value like ^ leads to an effect that can be estimated theoretically and measured 
experimentally. It was pointed out in reference [T7] that sizeable direct CP-violation in 
D° — > (p°,w)7 can be induced through Im[CV] provided that the LD amplitude carries 

2 Note this is the sign convention of [T5] but opposite to references [TBI ITT] . 
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a strong phase. The latter is necessary as the SD contribution of 0^ does not have a 
strong phase. Our work improves on [T7] in that it includes the (9g matrix element per se 
which carries a strong phase and thus does not rely on a size strong LD phase. The latter 
vanishes in the chiral limit at leading order for WA which we believe to be dominant. By 
how much corrections change the LD phase is an interesting open question for which we 



refer the reader at this stage to the conclusion and appendices A. 1.1 and B.3 



Other channels and effects that were proposed are the electric dipole moment of the 



nucleon [TBI US], CP-asymmetry in D° -»■ </> ->■ K+K~ [19j and D° V(-t PP) -»■ 
K + K~ [20]. We argue that the strong phase of the LD contribution, which we believe to 
be WA, is small due to chiral suppression at leading order in appendix [A] 

Another question is what type of NP models could induce such values as in ([5| without 
violating existing constraints. Among the NP models inducing (|5p are supersymmetric 
models [T31 ITB| l2~Tj . Randall-Sundrum flavour anarchy [22] and models of partial com- 
positeness [23], whereas in fourth family models it seems more difficult to accommodate 

EE]. 

The paper is organised as follows: In section [2] we introduce notation and remind 
the reader of why CP-violation in a charm system is supposed to be small in the SM. 
Section [3] is the main part of this paper and consists of detailing the two amplitudes 



3.1 and an estimate of direct and time dependent CP-violation in section 3J2 using the 
matrix elements of the operator [21] computed within light-cone sum rules (LCSR). 
Conclusion and discussion are presented in section |4j An important part of our work is 
the discussion around the estimate of the LD contribution. In appendix [A] it is argued 
that the WA contribution is the dominant LD effect and we comment on the possible size 
of the strong phase in A. 1.1 Appendix [B] includes a discussion of the SM CP-asymmetry 
in 



B.2 as well as the a discussion of the CV-effect, noted in [T7], in B.3 



2 Effective Hamiltonian and amplitudes 
2.1 I AC I = 1 Hamiltonitan 

Following, closely, the notation of [13] we write the effective AC = 1 SM Hamiltonian as 
follows 

T-L cG = XdHd + X s T-id + AfoHpeng , Ad = V* d V u d , D = d, s,b (6) 
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and 



U q = %\ ( "'O'! ■ h.c. . <! = <!.* 



* i=l 



x b n pcng = ^ (c 7 o 7 + C' r O' 7 + C 8 Os + c 8 o' 8 + ...) (7) 

with = 7 M (1 — 75) and a, (3 being colour indices. The Hamiltonian "H peng contains all 
the SD transitions including electric (A. 10) and chromomagnetic ^ operators as well as 



the four quark operators with structure different from 0\p- As compared to [15J we have 
absorbed the into the Wilson coefficient which is non-standard for the SM contribution, 
but convenient for our presentation as we can write: 

CP = CP SM + CP NP (8) 

with Cf p as in the literature. Since A djS = 0(A) and X b = C(A 5 ), where A ~ 0.226 pQ is 
the Wolfenstein parameter [T], one gets using the unitarity relation 

X d + X s + X b = 0, X d ~ -X s , X b ~0, (9) 

where the symbol ~ above is to be understood as up to corrections of 0(X 4 ). The fact 
that the third generation decouples up to 0(X A ) is the reason why in the SM the generic 
expectation for CP-violation is Aqp — few • 0(X 4 ) as mentioned in the introduction. Hier- 
archies in amplitudes might lead to an enhancement or further suppression. Throughout 
this paper we shall implicitly employ the Wolfenstein parametrisation. 
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2.2 D — > V^j transitions 

In this subsection we shall introduce some notation which aims to ease the presentation 
in the following subsection. We write the amplitude as follows [^J 

A[D ^ Vj ] = ( V j\n eS \D) = Ay + Ay , 

+ (1 0) 

witffl 

p i = 26 pa/37 6*V>V , P2 = 2i{{p- q )( v *-e*)-(r)*-q)(p-e*)} , (11) 

where r)(p) and e(q) stand for the vector meson and photon polarisation tensors respec- 
tively. The two structures 

A L{R) = (A 1 ±A 2 ), (12) 

of negative mass dimension, correspond to left- and right-handed polarised photons. The 
rate [33], in our conventions, is given by 

r[D -+ V>y] = ^-ml (l - 4V (IA| 2 + |A| 2 ) • (13) 
327r \ mj) J v ' 

3 CP-asymmetries in D — > V7 

The operators ^ consist of c — > u transitions of the FCNC type. In a heavy-to-light 
transition for which LCSR can make predictions j2l] the c-quark can pair with a u, d or s- 
quark leading to the transitions D° —y (p°, u) r j ) D + — y p + 7 and D£ —y K* +r y respectively. 
These transitions shall be investigated. Note it is only for the neutral D°-system that 
oscillations and thus TDCP asymmetries are feasible. 

As previously mentioned and outlined in appendix [C] direct CP- violation originates in 
its minimal form by two amplitudes with weak and strong phase difference. These two 
amplitudes will be given by the LD and Og contributions respectively, to be discussed 
below. 

3 The amplitudes A1.2 up to phases are often denoted by „4pc,pv m the literature e.g. [33J[T7]. The 
acronyms PC and PV stand for parity-conserving and -violating respectively. 

4 The sign convention for the epsilon tensor is given by tr[7 5 7 a 7;,7 c 7 ( 2] = 4ie a (, c d and are the ones used 
in the classic textbook of Bjorken & Drell. 
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3.1 Amplitudes with strong and weak phase difference 

LD: The LD contributions are generated by four quark operators amongst which are 
the WA (Figjljleft) and the quark loop (QL) (Fig|IJmiddle;right) topologies. We 
shall summarise the essential points below and refer the reader to appendix [A] for 
more details. The only global assumptions we need for our argument are that the 
LD contribution is a) dominant, b) has a small weak phase and c) a small strong 
phase. Points a) and b) are a direct consequence of the CKM/Cabibbo-hierarchy: 
WA terms are proportional to \d, s — 0{X) which has a small weak phase 0(X 4 ) and 
is large compared to the SD part which is proportional to A& ~ C(A 5 ) Q. 

Turning to the smallness of the strong phase let us summarise the essence of 
appendix [Aj Branching ratios are compatible with WA-dominance from LCSR- 
predictions of WA diagrams [26J. QCD estimates suggests that WA dominates 
QL by almost or about two orders of magnitude, as the latter are suppressed by 
two additional loops with respect to the former. We thus expect the lion's share 
of the strong phase to originate from radiative corrections. Further discussions are 



deferred to appendix A. 1.1 



We extract the amplitudes for D° — > (<j), K*°)~f from experimental branching ratios. 
The latter differ from LCSR predictions [26] by about a factor of 2 which translates 
into a factor y/2 ~ 1.4 on the level of the amplitudes which is rather good in view 
of the experimental and theoretical uncertainty. The latter is mostly due to radia- 
tive corrections, higher twist contributions and (outdated) distribution amplitude 
parameters. Scaling the LCSR-predictions of [26J by a factor a/2 we get with (A. 5) 
and ([12 



. . -5.8- 10" 8 . . -1.6 - 10~ 8 r ~ . 

A L ld ^ , An : ld , {D ->• 0°, u)} , 

a , 12.5 -10- 9 . . 0.7 -10- 9 rr ^ + 

A L ld ^ , A R ld ^ , {D+ -)• (p + ,K*+)}.(U) 

c v m D c v m D (a > b> 



5 As previously emphasized we do not distinguish between the p° and oj as the effect largely cancels in 
the ratio. 
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0$: The amplitude of the chromomagnetic operator (|3| is parametrised fcryj^J 



V27r 2 7 c y 1 (C 8 -QG 2 (0) 2 = 2 



where "H |s = -^(CsOg + C' 8 0' 8 ). Therefore Gi, 2 (0) corresponds to the matrix 
elements, with on-shell photon q 2 = 0, 



(Vl\0P\D) = (^) -(G!(0)Pi ±G 2 (0)P 2 ) , (16) 



which are analogous to the penguin matrix element for Ti and T 2 Eq. (A.ll). The 
variable e = y/Ana > is the electromagnetic charge. These matrix elements have 
been computed with LCSR in [22]. At leading twist-2 it was found that Gi(0) = 
G7 2 (0) which implies that 0$ and 0' 8 generate solely left- and right-handed ampli- 
tudes respectively. Moreover Gf°^ p ° 7 (0) ~ Gf°^ 7 (0), G? s+ ^* +7 (0) ~ Gf + ^ p+7 (0) 
to an accuracy sufficient for our purposes. We shall therefore not distinguish be- 
tween themQ The imaginary part, relevant for the CP-asymmetry, is found to be 

Im[G?°(0)] ~ -0.20(8) , Im[Gf + (0)] ~ -0.10(4) , (17) 



where numbers were rounded. The values in (|17j) are sizeable compared to typical 
estimates Tf° (0) ~ T ] D+ (0) ~ 0.7 of the-07 operator as compiled in [I7j . The 



difference in the numerical value of neutral and charged matrix element in Eq. (17) 
originate from the charges of the valence quarks of the mesons. Using the reference 
value for Im[Cg ] the relevant ratios | Al, 2 |sA<4i,2|ld| are around 10 -2 and thus the 
scale for CP-violation is set at the percent level. 



3.2 Direct CP-violation 

Since the photon polarisation is not easy to measure in practice a slightly inclusive rate 
T[D — > Vj] = T[D — > Vjl]+T[D — > Vjr] is measured. We parametrise the corresponding 

6 The factor cy is inserted to absorb trivial factors due to the uj ~ (uu + dd)/y/2,p° ~ (uu — dd)/y2 
wave functions, cy = — v2 for p° in c — > d, cy = y2 in all other transitions into ui & p° and cy = 1 
otherwise. Note in the overall CP-asymmetry this factor will drop out. 

7 In fact the ratio of the WA to the Gi(0) form factor is well approximated by R = r p /r u where 
rx = (fx)/( m xfx) i s the ratio of the tensor to the vector decay constant. Information on this ratio 
exists only sparsely in the literature. Similar remarks apply to the Dj —¥ K* + and D + — > ^-transitions. 



7 



amplitudes as follows, 



A l ,r = Ai±A 2 = h, R e i&L ' R + g L , R e iAL - R e 



t<PL,R, 



with 



lL(R) 



\h =t Z2I , h,2 = <A.i,2\ 



LD 



Gi(0) 



i A 



|(-) l |C «| (2Gl( 0)) 



ICJe 



\C' 8 \e i4 > R 



(19) 



where A i J? , and 4>l,r are the strong and weak phase of (15) respectively leaving the 
quantities Il(R),9l(r) real-valued. In the equation above we have made use of Gi(0) = 
£2(0), found at leading twist [M] . implying that 0% and 0' 8 solely contribute to the left- 
and right-handed amplitude respectively and in addition leads to A^ = A R . The latter 
is not true when the contribution due to ImfCy ] is included, in which case the formulae 



f° r 9lr have to be modified according to Eq. (A. 26) in appendix D 



In the case where the two photon polarisations are not distinguished the formula for 



CP- violation is slightly more complicated than the one given in Eq. (A.25). The general 



formulae and a derivation, including TDCP-asymmetries, can be found in Ref. [25] for 
example. We find 



A C p(D° Vj) 



-4 



n 



(g L l L sin(A L - 5 L ) sin(0 L ) + g R l R sin(A R - S R ) sin((p R )) 



n = 2(ll + l 2 R + 2 (g L l L cos(A L - 8 L ) cos(0 L ) + g R l R cos(A jR - 8 R ) cos (</>#)) + g\ + c 

Assuming Il{r) 3> gL(R) and imposing A = Al = A R one gets: 

-2 



Aqp(D Vj) 



, , ?2 {gdh sin (A - 5 L ) sin(</> L ) + g R l R sin(A - 5 R ) sm((j) R )) .(21) 
l l + l R 



In the absence of a computation, and in view of the chiral suppression at leading order, 



we shall set the LD phases 5l >r (18) to zero in remaining formulae. This allows us to 
express A CP in terms of quantities discussed at the beginning of the paper: 



*Vy) = ^ (5) ™ ( Wm[C8 , + «]) . 



(22) 
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This formula, modulo notation, reduces to Aqp (A. 25) for l\ = l 2 (i.e. Ir = 0). 
With m c = 1.3 GeV, Eqs. (J5]), (14) and (17) we get for the neutral transitions 



A CP (D° -> (p°,u)i) = (-3.84Im[Cf p ] - 1.0Alm[C' 8 NP ]) c B 

( Im[C 8 wp 



-1.5% 



o.4% (; m i c " p i 



.4-10-V VO.4-10- 2 
with an estimated of uncertainty of about 45%, to be discussed below, and 



(23) 



cb 



cb 



0.6 x 10- 



1/2 



(24) 



B(D° -> (p°,w) 7 ) 

being the correction factor for the yet to be measured branching ratios. In going from 



(21) to (23) we have used the fact that the imaginary part of Cf M , which contains the 



CKM prefactors, is negligible with respect to the values (J5|. For the charged transitions 
we get 

^cpP+ s) (p + ,iT + ) 7 ) = (9.71Im[Cf p ] + 0.60Im[C^ p ]) c B 



3.9% 



( lm [Q 



NP] 



0.4 • 10 



-2 



0.2% 



V0.4- 10 



NPl 



(25) 



with an estimated of uncertainty of about 45% to be discussed below. Note, the dominance 
of the CVcontribution, in both neutral and charged case, is due to the Al(r) hierarchy in 



Eq.pi]). In fact the different sensitivity of AA C p, A C p(D° -)> (p ,^) and A C p(D 



(d,s) 



(p + , K* + )'y with respect to Im[Cs] and Im[Cg] gives a handle to discriminate between the 
individual contributions of the two chromomagnetic operators. 

Let us turn to the discussion of the uncertainty. The major uncertainty comes from the 
estimate of the 0& matrix elements which we estimate to be around 35% |24j . Then there 
is the phase of the WA contribution, 5l,r, for which we assign an uncertainty \Sl,r\ = 10° 



c.f. appendix A. 1.1) which leads to a 20% uncertainty. Amongst the LD contributions 



the combination l\ + l R is taken from experiment but the ratio Il/Ir which we took from 
pZB] could have uncertainties, say, at the 20%-level. Adding the three sources discussed 
above in quadrature, as they would seem uncorrelated, we get about 45% uncertainty. A 



few additional remarks are in order. In appendix |B.2| we estimate the SM contribution to 
be of the order of 10 -4 which is negligible. Furthermore we refrain from including at this 



point the uncertainty due to the C 7 -effect discussed in appendix B.3 We would like to 



9 



mention though that it cannot be excluded, depending on the model and the LD phase, 
that the C-j and Cg-effect conspire to cancel significantly in the CP-asymmetry. 



3.3 Time-dependent CP-violation 

As a result of D°-D° oscillations CP-asymmetries are time dependent for the neutral me- 
son, which will lead to novel features. In particular TDCP asymmetries do not necessitate 
a strong phase difference in the two amplitudes. Thus in principle we have to adjust the 
amplitudes to include the CV-effect, from Eqs. ( l~8fl~9 ) to ()l,r (A. 26) as detailed in 
appendix [D] Indications are though that these effects are overshadowed by the dominance 
of the LD amplitudes Il,r- 

Important mixing parameters of the D°-D° system are the mass and width difference, 
the mixing phase 4>d as well as the ratio \p/q\ of the parameters p and q translating 
between the flavour and mass eigenstates. The latest HFAG values [27J are 



xd 



Amp 



0.62(19) • 10" 



=1.04(7) 



D 



Ar 



Vd 



D 



2T 



0.75(12) ■ 10" 



!>d = - 10.1(9.0)°[-2.02(2.70) c 



(26) 



where V = (r^o) -1 is the inverse lifetime of the Z)°-mesons and Attiq and Ar^ are the 
difference of the heavy and the light D°-meson mass and width respectively. We quote 
the values with no direct CP-violation except for (pp, where we also quote the result 
for allowed direct CP-violation as it is rather different in the central value, although 
compatible within errors. 

In view of the closeness of \p/q\D to unity in (26) we shall assume \pjq\n = 1 in which 
case the TDCP-asymmetry assumes the following form, 



A CP (D^ Vi)[t\ 



S sin( Am ot) — C cos(Ami)t) 



cosh(^t) - Hsmh(^t) 



(27) 



where the convention Aqp{0) = —C is somewhat awkward but standard. The formulae 
for S and H are given in appendix [D] and C = — Acp from the previous section. Let us 
define the LD chirality-asymmetry (ratio) by 



Xld 



/ 2 



'"2 



21 T I 



if + q 



il + il 



(2f 



Note for the values in (14) we get xld — 0.5(1). Thus if we assume xld ^ 10" 



10 



Il,r 3> 9l,r, which both seem true, and once more set 5l,r = we get an interesting 
expression for for H and S, 

H W-W7% -(-ecos[sin](0 D )) 

l L + L R 

= XLD-(-ecos[sin](0 D )) , (29) 

which directly measures the ratio of the LD chirality structure times the cosine and sine 
of the mixing angle of the D°-system. The variable £ = ±1 is the CP-eigenvalue of the 
y-meson whose values can be found in appendix [Dj With £(p°,u) = 1 we get 



H[D° -> (p°, w) 7 ] ~ -0.5(1) cos(0 D ) , S[D° ->■ (p°, w) 7 ] ~ -0.5(1) sin(0 D ) . (30) 

Let us emphasize once more that this relation is valid in the case where a left- and right- 
handed amplitude are of comparable size and dominate all the other contributions. These 
ratios do not depend on the branching ratio and therefore do not have a correction factor 



c B (24). 



The experimental tractability of S and or H depends on the angle </>£>. Should (pc (26 ) 



that is to say sin0£>, turn out to be sizeable then S could be measured as for B — > K*^y 



at the B-factories. If coscpo is sizeable, which is what the value in (26) indicates, then 
one would need to focus on H. The latter might be measured, in analogy to B s — > (jyy 
case [25J, in the rates D° — > (p°,oj)j and the one for D° without flavour-tagging, which 
has experimental advantages, though it has to be added that the relatively small width 
difference in the D° system, yD/yB s — 0.1, means that roughly a hundred times more data 
has to be accumulated to achieve the same precision on H in the D°- as in the -B s -system. 

4 Conclusions 

Partly building up on ideas in [17] we have shown how Im[Cg ] become observable in 



CP-asymmetries in D — > Vj. Setting the LD phases 5l,r = (18), in the absence of a 
computation, we got 

A cr <D{ d „ } -» A- h ) . (3.9% gffi) + 0.2% (0^) ) c B , (31) 
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with eg (24) a correction factor for a yet to be measured branching fraction of D° — > p 0/ ~f. 



Uncertainties are in the 45%-range, c.f. section 3.2, for which we have not included the 
uncertainty of the Im[CV]-effect on which we comment in the the last paragraph of this 
sect ion. The SM contribution is negligible, down by two orders of magnitude c.f. appendix 



B.2 



for reference values of Im[Cg ^ NP ] used throughout this paper. An obviously interesting 
aspect is that the Wilson coefficients of the two chiralities of the chromomagnetic operator 
enter with different sensitivity. Even the relative sign is different between the D — > V'-f 
and the D° -)• KK, ttit case (Ji]), 

AA™p ~ -1.8 (lm[Cf p ] - lm[C' 8 NP }^ sin(5)c NF , (32) 

where cnf is a correction factor for using naive factorisation of order one. Note whereas 
the sign of the entire contribution is predicted for D — > V'y this is not the case for the 
non-leptonic mode as the sign of sin(<5) remains undetermined. 

The chirality of the photon is an interesting aspect and deserves some discussion 
in comparing it to the 6-sector. In b — > (d, 5)7 transitions the left-handed amplitude 
dominates over the right-handed amplitude as a result of the large 6-quark mass and 
the V-A interactions. This pattern might be broken by physics beyond the SM and can 
be measured in TDCP-asymmerties [28]. The situation in D° — » V7 is rather different. 
Whereas it is still true that the left-handed amplitude is larger than the right-handed 



amplitude, e.g. (14) it is not very significant since the c-quark mass is smaller. This 
neither-nor situation has two consequences. 

First since the amplitudes themselves are LD dominated the TDCP-asymmetries are 
not sensitive to novel right-handed currents. On a positive note TDCP-asymmetries mea- 



sure the LD chirality asymmetry xld (28) and thus can provide interesting information 
on LD dynamics and could serve as validation criteria for theoretical tools. Let us add 
that the feasibility of the measurement depends on the definite value of the mixing phase 



d as commented on at the end section 3.3 Second the fact that the CP-asymmetries in 



D — > V'y are more sensitive to C% than C' 8 (31) follows from Al > Ar, c.f. Eqs. ( 14|12 ). 



Thus the direct CP-asymmetry in D — > is not a good place either to look for right- 



handed currents as reflected in (31) in the low relative sensitivity of C 8 versus C\ 



On the speculative side it is of course possible that NP contributes to SM or non-SM 
operators of the WA-type, Of 2 



possibly with new weak phases. Allowing for the 
latter and parametrising a strong phase for the yet to be computed (9(a s )-corrections 



8 In addition, extensions of the SM which remotely follow the pattern of minimal flavour violation 
|Cg/Cg| ~ m u /m c predict small effects in C 8 per se. 

9 Note in [15] it is the GIM-combination ([9]), Of 2 ~ Of 2 i which is severely constrained through e'/e in 
new weak phases but not the individual operators °f down and strange per se. 
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k(R) -» ii(ij)e i$i(fl) , one gets: 

= Xld ■ (-£cos[sin](0 D - $ L - cos(5 L - (33) 

Needless to say that xld is then affected by the NP. 

At last let us give an outlook and hint how the current work could be improved. On the 
theoretical side it would benefit from a thorough reassessment of the WA contributions, 
that is to say the full (9(a s )-correction and using updated values of the photon DA [29] 
and other parameters. In particular the radiative corrections would allow an estimate of 



the strong phase and the inclusion of the CV-effect jTT], c.f. appendix B.3, to the central 



value rather than as an error estimate. Our discussion in appendix A. 1.1 suggests that 
the convergence of the a s -expansion is at least good enough to get the sign and a rough 
idea of the value of the angle. This would allow exclusion of scenarios where the C-j- 
and CVeffects conspire to cancel each other in the direct CP-asymmetry, c.f. appendix 



B.3 The prominence of WA in the isospin asymmetry in b — > s processes provides yet 
another motivation for their reassessment. Furthermore it might be interesting to extend 
this work from D — > V^f to D — > Vl + l~ as the latter might be easier to deal with at 
the LHCb where the photon final state remains challenging at present. The estimated 
branching ratios for D° — > (p°,u) , ~f are a factor of 30 below the current experimental 



limit 10 Should the CP- violation in non-leptonic charm decays ([TJ remain one would hope 
that A CP (D — > Vj) is going to be measured at some future flavour factory or possibly at 
the upgraded LHCb. 
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A Long-distance contributions 



In section 



A.l 



we extract the LD contributions in D° — > (p°,u;)7 and Dt d , — > (p + , K* + ) , y 



from experiment with minimal theoretical input. In section A. 2 we comment on the 



10 



For the D^ d ^ — > K* + )j decays there are not even limits, possibly because they are experimen- 
tally more challenging. 
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hierarchy between the weak annihilation (WA) and the quark loop (QL) contributions. 



A.l Estimate of LD contribution 



Our basic assumption is that WA, as discussed in section |AT2| dominates over QL. There- 
fore we feel justified to use LCSR result in as a guideline for the ratio of left- to 
right-handed amplitude as well as the sign of each of them. More precisely we infer that 



the strong phases are relatively small, to be discussed in more detail in section |A.1.1 
and thus we assume the central value to be real valued. This information is sufficient to 
extract the amplitudes from experimental branching ratios. 
The measured rates of the D° — > V'j-tjpe are [TJ, 

B(D° -)• iT 7 ) = 3.27(34) • 1(T 4 , B(D° -> <fyy) = 2.70(35) • 10" 5 , (A.l) 



of the Cabibbo allowed and singly Cabibbo suppressed type. From the rate (13), with 
t d o ~ 4.1 • 10~ 13 s, we get 

x[K*] cx = 1.8 • 1(T 14 (10%) GeV~ 2 , x[0] ex = 2.0 ■ 10~ 15 (10%) GeV~ 2 (A.2) 

where we have introduced the following shorthand, x[V] = (|Al| 2 + \A2 | 2 ) £) o_ s .y 7 , which 



is the rate (13) modulo phase space factors. The relation between Aip POj) and A P c(pv) 
of reference [26], for cu — > dd (e.g. D° — > (p°,w)7) transition, is given bwH 

Ai(2) |wa,o(o°) — — f^^ dai ) \ ~c V ) ^ PC(PV ) (^--3) 

where a x — C\ + C 2 /3 ~ —0.5 is the colour suppressed Wilson coefficient. The ratio of 
the LCSR prediction to the experimental value is 

Cj^ - °- 5 < A - 4 » 

x[K*\ ex 

which translates into a factor of v2 for the amplitudes. Thus the LCSR predictions 
of WA saturate the central value of the experimental result up to 40%. In view of the 
experimental and theoretical uncertainty this is a rather good result. The theoretical 
uncertainty is due to hadronic input parameters, omission of radiative corrections and 
semi-global quark-hadron duality. We shall therefore take the predictions of for 



11 With respect to reference [26\ the labels 1 and 2 are chosen the other way around in % eff Q. 
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D° -»■ (p°,w)7 and -> (p + ,-^* + )7 and scale them by a factor of — l.^Pj In 

doing so we implicitly rely on the sign and the ratio both of which would seem more 
reliable than the individual amplitudes. In following our recipe above we get, 

A(2)|lD - 1-4 X -4l(2)|wA,0(aO) , (A.5) 

using the notation = Ai(2) |ld we get 

, -3.7 -10- 8 -2.1 - 10- 8 rn0 

Cym D CyVTlD 

In turn this gives a result for the branching ratios of 

£(D ^(p°,u;)7)| S ~0.6-10- 5 , 

£?(£>+ iT + 7 ) ~ 2i3(D+ p + 7 )|jX5l - 0.5 ■ 10~ 6 (A.7) 



Comparing — )■ (p°,o;)7) above with B(D° — > 07)pdg (A.l) we note that the former 
is down by about a factor of four. A factor of two comes from Cy = 2, i.e. the valence 
quark content of p° and u. Another factor 1.2 comes from phase space factors and myfv- 
prefactors. The remaining discrepancy might partly be due to neglecting the strange 
quark mass. It might be questioned whether 2m s , which results from the axial current, is 
really that small a parameter. At last we note that the value for the p°-branching ratio 
is well below the current limit of 2.4 • 10~ 4 pp. 

A. 1.1 Strong phase in WA-process 



Here we shall add comments on the the strong phase Slr (18) in WA, that is to say, 
we discuss the suppression at leading order and reflect on the reliability of perturbation 
theory for this computation. 

The photon can either be emitted from the initial state meson or the final state meson, 
c.f. Figjljleft) . At leading order 0(at®) the former does not lead to an imaginary part 
and the latter is chirally suppressed. More precisely the vector part vanishes identically 

12 We should add that this itself does not settle the issue of the hierarchy between the WA and QL, to 
be discussed in subsection |A.2| ,as the decays above do not have QL or at least if the uu component of 
the 4> is neglected. 
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by vector current conservation and the axial part, which corresponds to Ai, is partially 
conserved; that is to say proportional to 2m ( i and thus negligible as noted in [2SJ Eof^j 

The reason we can trust perturbation theory, beyond questions of convergence of the 
series, is that m 2 = and m 2 D ~ (1.86 GeV) 2 are sufficiently far away from the narrow 
peaks such as the u (m u — 0.78 GeV). One can get an idea of this by looking at the 
famous cr(e + e~ — >■ e + e~) cross section plot [I] in the region of y/s ~ for example. 
Note the latter probes the electromagnetic current which is not the same but definitely has 
similarities to the ^7^(1 — 75 )d current relevant for WA. The partonic over the hadronic 
cross section at mo indicates about a 25% correction and the 3-loop perturbative QCD 
result is roughly in between the two. Thus it seems fair to say that in the region of 
perturbative QCD is neither great nor bad, somewhere in between. Let us now turn to the 
question of the size of the (9(a s )-corrections from another point of view. The (9g matrix 
element might be regarded as a radiative correction to the O7 matrix element, at least the 
part where the gluon does not connect to the spectator quark. The latter contribution is 
about 1/2 and 2/3 in the neutral and charged case respectively [21]. Taking the average 



of the charged and neutral (17) matrix element multiplying by the relevant factor from 
above one gets (0.2 x 1/2 + 0.1 x 2/3)/2 ~ 0.085 which has to be compared to the reference 
value Ti(0) = 0.7 chosen throughout this paper. This would indicate C(a s )-correction 
around the 15%-levef^l 

In view of the discussion above we shall take 25% as a, probably conservative, reference 
value for radiative corrections. If we assume that real and imaginary part share out this 
value in an equal way we end up with \5l\ — 10° and \8r\ ~ 10°[^} Needless to say 
that a computation of the radiative correction would add more weight to these kind of 
discussions. One might be sceptical about the convergence of the a s -expansion in the 
mo region yet it would seem that the first correction should at least give a reasonable 
estimate of the phases 5l ft and in particular of their directions which is relevant for 



excluding cancellation effects between the C 7 - and C 8 -effect as discussed in appendix B.3 



13 Note there are also four quark operators other than O1.2 whose matrix elements are not chirally 
suppressed but since they originate from SD physics their Wilson coefficients are tiny in the SM as 
previously discussed. 

14 There are two further aspects that one could bring in. Only a limited number of graphs generate an 
imaginary part at (D(a s ). On the other hand the radiative correction are colour-enhanced as with respect 
to tree graphs. Possibly these two effects just balance each other off so we shall not give it any further 
thought here. 

15 The two phases ought to have the same sign as a remnant of the heavy quark chirality structure. 
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A. 2 Weak annihilation versus quark loops - theory point of view 



One may distinguish two types of LD contributions according to the quark level transitions 
is cu — y dd or c — y udd(ss) generated by the weak operators Of'^ for instance. From 
the viewpoint of quarks and gluons the first type is known as weak annihilation (WA) 
(Fig. [TJ left) and we shall name the second type quark loops (QL) (Fig. [TJ middle;right). 
The WA contributions have been computed in 1995 for B —y V'j and D —y Vj in [3T] 
and B° — >■ (p ,oj)j in [2H ED] at O(a s ). We hasten to add, for further understanding of 
our arguments below, that QL of the type shown in Fig. []J middle, right) are evaluated in 
an l/m c (l/mfe)-expansion for c(b) —> u(d,s)j, although in principle one could compute 
them in the exclusive case with LCSR which does not adhere to a l/m c (l/m&)-picture. 
We shall argue that we believe the WA to dominate over QL: 

• Generic observation: QL and WA are generated by the same weak operator, Of'2 
and Of 2 @ respectively, yet the QL is down by two loops with respect to WApj 
This is the case because the single quark loop Fig. [l]( middle) vanishes exactly by 
virtue of gauge invariance. The photon polarisation H^q) = {q 2 g^ u — <Z/i<Zw)n(g 2 ) 
vanishes for q 2 = when contracted with the photon polarisation tensor. This 
suggests a natural hierarchy WA ^> QL in the types of charm transitions discussed 
in this paper. 

• Test-case in B-physics: One would think that the perturbative picture is valid in B- 
physics. Taking numbers from [32] for the WAp^jand QL one gets: \Aql/Awa\b 
2 ■ 10~ 2 . To be more precise, for Aql we have taken the charm loop contribution 
where the gluon is radiated into the final state vector-meson. The latter, as men- 
tioned above, does not depend on the l/m^-expansion. Note WA for B° —y p 0/ y is 
accidentally small because of cancellations between tree-level and penguin four quark 
operator contributions. We do not expect the same to take place for D° — y (p°,u;)7 
since those cancellation are between tree and penguin four quark operator contri- 
butions and the latter are tiny in D-physics. 

• Test-case in D-physics: Does this hierarchy remain intact for the D-physics. De- 
spite the obvious fact that the a s (m c )-expansion and the l/m c -expansion are less 
trustworthy it seems hard to see how a hierarchy of two order of magnitude can be 
overthrown. Taking the contribution Fig. [T](right) for the QL from |31j . which does 



p 7 



16 In principle there is a GIM suppression of the QL in addition which is though not very effective for 
the matrix elements [5T] . 

17 WA is Cabibbo suppressed with respect to QL in _B-physics. In comparing the WA and QL pro- 
cesses/diagrams we, of course, do not take CKM hierarchies into account, especially because they are not 
present in the charm decays we are interested in. 
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rely on l/m c -expansion, and the estimates of [26J one gets a number, \Aql/Awa\ — 
2 • 10~ 2 , which is surely accidentally close to the one for the B~ — > p^'y. An estimate 
of the size of l/m c -corrections, for WA, can be gained by looking at the chirality 
structure of the D(B) — > Vj result in LCSR. In the heavy quark limit it is believed 
that Ai = A2 (i.e. Ar = 0). The results in [26] indicates that this is fulfilled for 
the D°(B°) — > p°7 transition at 5 7% (70%) -level respectively. This does not suggest 
a dramatic breakdown, by which we mean one or two orders of magnitude, of the 
1 / m c -expansion. 

Thus our analysis suggests that WA dominates QL by roughly two orders of magnitude. 
We shall briefly comment on another approach. In the extensive work [33] the two transi- 
tions were modelled with hadronic data. It would seem that WA corresponds to the pole 
(P) terms and QL to the vector-meson dominance (VMD) part. Comparable numbers 
for P and VMD were found which is not in line with the arguments abovd^J One might 
wonder how the vanishing of the single quark loop in Fig. [T]( middle) is reflected in this for- 
malism. A problem is that the signs of the couplings of the VMD models are not known, 
that is to say their absolute values can sometimes be taken from experiment. Thus the 
formalism might overestimate the contributions as it cannot capture cancellations, which 
gauge invariance almost suggests to be present. A similar point of view has been taken 
in [31] by one of the authors of [33 J in chapter 3.1.3. 




Figure 1: A selection of LD diagrams for D° — > (p ,uj)j. Note it is the fact that the p°/u carry 
both dd and uu components that makes the same operator 0^ s ([7]) contribute to both (WA & 
QL) topologies, (left) Weak annihilation (WA). (middle) Quark loop (QL). This contribution 
vanishes, exactly, for on-shell photon by virtue of gauge invariance as discussed in the text, 
(right) QL example of C(a s )-correction. This diagram has a sizeable imaginary part which can 
be inferred from the computation for c — > wy in reference |3T) . 



18 We note that in [53] the P-part receives no contribution in .42 (o Apy) which is a fact that is not 
reflected in the LCSR computation [251 [50] , 



18 



B A CP (D° -> I/7) other than through C^ p 

For our discussion it is convenient to write the amplitude as follows, 

A ~ A d e i5 M d + A a e w M s + X b e i5b A b , (A.8) 



which is similar to (A.24) with the exception that the unitarity relation ^ has not been 



used and that the weak phases are contained within Xd, s ,b (j6])- As argued in appendix |A| 
we expect the lion's share of to be covered by WA which has, presumably, a small 
strong phase which we shall neglect (Sd —> 0). We assume a Wolfenstein parametrisation 
up to order C(A 5 ) which fulfils, e.g. [35] . 

Im[A d ] = , Im[A s ] = A 2 X 5 i] , Im[A 5 ] = -A 2 X 5 i] (A.9) 

where A, p and 77 are the other three Wolfenstein parameters and A 2 X 5 r] ~ 1.4 ■ 10 . 



Eq.(A.9). The fact that |Im[A& )S ]| ~ 1.4 • 10~ 4 indicates small CP-asymmetrie^j, of that 
order. 

Thus it remains to identify contributions with sizeable strong phases 5 St b and ampli- 



tudes A s fi for which we see two major sources. First the matrix element of Og e.g. (17) 
and second the matrix element of 0^ s [31] (c.f. Fig. [l]( right) for a contribution) giving 



rise, effectively, to an CVoperator. The latter as well as its matrix element analogous to 



(16) are defined and parameterised respectively as follows, 



{ 7 ] = -^ua^(l± l5 )c (A.10) 
<Y>y\oP\D) = (^) ^(T 1 (0)P 1 ±T 2 (0)P 2 ) . (A.ll) 
B.l Effective Wilson coefficients Cj S 8 (m c ) 

Let us state that we do not intend to give a critical review of the treatment of Wilson 
coefficients in the charm sector, e.g. of whether it makes sense to include light-quarks 



into SD contributions evaluated in perturbation theory 20 We shall simply follow the 
literature. It is fortunate that the SD contributions turn out to be subdominant in the 
SM. 



19 One might be tempted to say that if WA dominates by another two order of magnitudes then this 
implies that the CP asymmetry is automatically below 10 -5 . This is not correct as in this way of thinking 
the absolute value of should be factored into Ab and then Im[Af,/|Af,|] ~ 0(1) is not small any more. 

20 We are grateful to Ikaros Bigi and Ayan Paul to draw our attention to this point. 
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The different contributions discussed above are conveniently discussed in terms of 
so-called effective Wilson coefficients. The latter consists of the pure Wilson coefficient 
CV,8( m c) an d matrix elements which can be rewritten in terms of O7 8 which we denote 
by 5C^ s 8 (m c ) 

C 7 cff 8 (m c ) = C 7 , 8 (m c ) + 5Cf 8 (m c ) . (A.12) 

From a conceptual point of view the Wilson coefficient can be divided into two further 
sub-parts, 

C 7i8 (m c ) = CgrVc) + C\j\m c ) . (A.13) 
The notation above is non-standard but hopefully useful for clarity. For the reminder of 

/ x 14 14 

this section we closely follow the notation of [36] . For C| ff only Cg (m c ) = r]c 5 r] b 23 C 8 (mw) , 
r/b = a s {mw) / C( s {mb) and r\ c = a s (mb)/a s (m c ), is known explicitly in the literature. For 

'efi 



C| ff all three parts are known which we shall quote, almost explicitly, below 



16 16 \Q / 14 14 16 16\ 



C\ mb \m c ) = -X b ^ C 3 [m h )X 3i ^ , (A.14) 



1 j 



where i = 1..8, j = 1..6. Note Cj mw \m c ) describes the evolution directly from my/ to m c 
and Cj mb \m c ) originates from integrating out the 6-quark at the m^-scale and running 
from rrib to m c . We hasten to add that the above expressions are given in the leading 
logarithm approximation. The term from the four quark matrix element is given by [31] 



5Cf(m c ) = ^Ac 2 (m c ) (\ s f[(m s /m c ) 2 ] + X d f[(m d /m c ) 2 ]) . (A.15) 

The strong phase results from the charmed meson's four momentum cutting the diagram 
through light quark lines. The contribution of Ci(m c ) vanishes whereas the C^^^(m c ) 
have not been given but are small as they originate from SD contributions which them- 
selves are small. In fact the numerical hierarchy is as follows [31]: 

\c\ mw \m c )\ ~ 2 ■ 1(T 7 < \C^ mb \m c )\ ~ 8 ■ 1(T 6 < \5Cf(m c )\ = 5 ■ 10" 3 . (A. 16) 

The hierarchy between the first two was noted in [33j and numerically improved in [31]. 
The fact that matrix element dominates the Wilson coefficient was pointed out in [31]. 
The expression of Cy" 1 ^ (jn c ) f° r operators other than 2 was given recently in ref. [32] ■ As 
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mentioned previously we are not aware of explicit results for Cg m ^(m c ) and 5C| ff (m c ) in 
the literature, yet they can be expected to be close to their Crcounterparts as they differ 
only by colour factors. Excluding cancellation effects we would expect them to equal up 
to 0(1/N C ) effects, say equal to about 30 — 50%. Given the uncertainties of the estimates 
the approximations, Cg (m c ) ~ C£ (m c ) and 5Cf ff (m c ) ~ 5C| ff (m c ), are good for our 
purposes^ Furthermore with C 8 (m c ) ~ Cg mb) (m c ) c\ mb \m c ) ~ (-0.3 + 0.8?) • 10~ 5 
we see that the SM value is two to three order of magnitude below the reference value 
Im[Cf p ] ~ 0.4 • fO- 2 . 

B.2 A CP (D° -> Vi) in the SM 

In the SM we identify three main sources contributing to the direct CP-asymmetry: a) 
Cg(m c ) ~ Cg\m c ) b) 5Cj S (m c ) and c) £Cg ff (m c ). Right-handed operators 0\\ are 
negligible in the SM as Wilson coefficients as well as matrix elements are suppressed. As 
previously mentioned we shall use Cs(m c ) ~ Cf{m^} for ) and c) which is good up 

to 1/N C corrections. Note, as the leading LD amplitude is proportional to Ad, it is only 
X s or A& that can contribute to the direct CP-asymmetry. 



a) It is found that [31] 

C\ mb \m c ) ~ 0.06A 6 ~ (0.3 - 0.8z) • 10" 5 (A.17) 

and assuming, as discussed above, Cg mb \m c ) « Cj mb \m c ), we get that this contri- 
bution compares with Cg in Aq? as follows: 

0.06Im[A 6 ] ^ n o in _ 2 



Im[C 8 NP ] 



0.2 • 10~ 2 . (A.18) 



b) It is found that 



SCf{m c ) = (0.6 + 2.2z) • 10~ 2 X s + c\ d , (A.19) 

where the imaginary part, other than A s , corresponds to a strong phase. The number 
c is of no importance for CP-violation as it can be absorbed into WA which is 
proportional to A^ and much larger. The contribution Aqp compares with Cg p as 
follows: 

im'\ limfin a ■■.■o oi) . in--l7\(o) ^ ^_ 2 (A 20) 



Im[C 8 NP ]Im[Gi(0)] 



21 Though the values Cy™ w '(m c ) differ substantially for various reasons but this is of no concern as 
they are small. 
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for reference values (j5]), Xi(0) = 0.7 and Im[G] 



-0.2. 



c) As discussed above we expect SC, 



eft. 



5Cj (m c ) and this leads to a result for 



c) with Im[(xi(0)]/Xi(0) — 2/7 suppression factor as compared to (A. 20). 



Summa summarum the SM contributions are two orders of magnitude below the values 
Im[C| JP ] (^-contribution and with the value in (23) we get 



A C p\sm(D° -> (p°,u)i) « (-1.5%c B )(-2 ■ 10- 2 ) « 3 • 10 



-4 



(A.21) 



with eg a correction factor for the branching ratio (24) which we set to unity in the last 



step. We refrain from quoting a specific uncertainty. We would though think that the 
value catches the right order of magnitude. As possible criticisms one could advocate for 



example the estimate C i 



( m b) , 



a 



( m b) . 



in (A.19) 



m c ) and question the accuracy of local duality 

The charged 



Im[Gf + ] in QA.20D and this would lead to 



We refer the reader to appendix [AXjj for related discussions 
case is obtained by replacing ImfGf 
3.9%c B (-3 ■ 10~ 2 ) : 



A \ D+ 
^cpIsm 



-1 ■ 10" 



B.3 Acp(D u — > V7) through Im[C^ p ] and a strong LD-phase 



In reference [UJ the idea was put forward that C 8 (???np) mixes into CV(ra c ), e.g. Eq. (A. 14) 
for the SM evolution. More precisely depending on the model and the scale of NP, Mnp, 
it was put forward [17J that this leads to comparable value^J An important point is 
that C-j{m^) hardly affects D° — > tttt/KK because of a-suppression and is therefore not 
constrained by the latter. Following [17] we shall assume only SM degrees of freedom 
below the scale Mnp = 1 TeV and that the NP part of the Wilson coefficients is much 



larger than the SM part. Amending the notation of (A. 14) to include the running of six 
quarks above the top threshold one gets 



a 



(1 TeV) 



m c ) w 0.42C 8 (lTeV) 



(A.22) 



Cf TcV) (m c ) w 0.37C 7 (lTeV) - 0.26C 8 (lTeV) w 0.37C 7 (lTeV) - 0.62C 8 (m c ) 



and the analogous equations for the 7 8 -operators. Eq. (A.22) exposes the dependence 
of CV(m c ) on the scale Mnp and C) (Mnp). We shall somewhat arbitrarily choose the 
value Im[C}' )NP (m c )] w -0.5Im[CP NP (m c )] as a reference values. This follows the model 
dependent assumption llmfcf^l TeV)]| < |Im[Cp(l TeV)] | in [UJ. 



22 Note our normalisation of C 7 differs from [17] by a factor of Q u which translates in to Q u Cj 
where IK stands for the authors of 1171. 



IK 
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Since the O7 matrix element itself, as opposed to 5C% S , does not carry a strong phase 
and the LD strong phase vanishes at leading order in the chiral limit, as discussed in 
appendix A. 1.1 we did not include this effect in our results ( 23p5 ). In fact we estimated 
that the phases could be around |<5l,r| — 10° and we shall investigate how the CP- 



asymmetry changes. It is then useful to rewrite the gi amplitude as in (A. 26) with the 
replacement: 

[C 8 (2G!(0)) + C 7 (2Ti(0))] £P pMCsKCMO) " 0.5^(0))] (A.23) 



For Ti(0) = 0.7 and Gf°(0) ~ -0.2 - 0.2z ~ 0.3e~ il35 ° [23] one gets F 1 ~ -0.55 - 0.2z = 
0.7e~ ll6 ° . Thus a correction of the LD phase 5l,r = ±10° leads to a strong phase 
difference between the two amplitudes in the range of 10° to 30° which corresponds to a 
rescaling of the CP asymmetry by factors sin(10°) / sin(20°) ~ 0.5 and sin(30°) / sin(20°) ~ 
1.5 respectively. Thus in conclusion one cannot exclude the possibility that the phases 
conspire to cancel a significant part, or even an order of magnitude, of the effect! A lot 
of things have to go wrong for this to happen though. As discussed in section | A . 1 . 1 an 
0(a s ) computation would presumably give an indication of the sign of the LD phase as 
well as its size and would allow to make firmer statements. 



C Formulae for direct CP-violation 

In this appendix we collect some formulae which are useful throughout the text. We shall 
parametrise an amplitude as follows, 

A(D° ->■ /) = A a e iSa e i<f)a + A b e iSb e^ b , (A.24) 

with weak (CP-odd) phases and strong (CP-even) phases 5 separated to leave A a>h real. 
Note that in the SM the decomposition (A.24) is sufficient as one might use unitarity ^ 



to eliminate one amplitude to arrive at two amplitudes. Using the notation A = ^ , 
8((j))ab = &{4>)a — 8(4>)b the CP-asymmetry becomes: 



Acp[D° -> /] 



-2 sin(5 a6 ) sin(0 ab ) 



1 + 2 A cos((U) cos(0 ab ) + A 2 

A<1 



-2sin(cU)sin(c/> ab )A . (A.25) 
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In the second line we have assumed a hierarchy between the amplitudes which is the case 
for D° — > (p°,u;)7 as studied in this paper. 

D Formulae for TDCP-violation 



The replacement due to the relevance of O7 as described in subsection 3.3 is as follows: 



9Le « e i*L g Le ^j* L = <?* 1^(2^(0)) + C r (2Ti(0))] , (A.26) 

and for g R is given by the following replacements: L — > R and C 8 ,C 7 — >■ C' 8 ,C 7 . Note 
unlike before we cannot assume a common strong phase as the ratios Cg/CV and C' 8 /C' 7 
might not necessarily be the same. This is why the strong phase A carries a chirality 



label. The symbol $ denotes the weak phase. The formulae for H and S in (27) are 
given, including a derivation, in the appendix of reference [25j^] and take the following 
form: 

H[S] = -^-[IlIrCOs(5 l - 5 R ) cos[sin](0 D ) + l L G R cos(A R - S L ) cos[sin](0 D -$ fl ) 

+ gdR cos(A L - S R ) cos[sin](0 D -$ L ) + g L g R cos(A L - A R ) cos[sin](0 D -$ L -$ jR ) 

N = 2 {l\ + l\ + 2l L g L cos(A L - 5 L ) cos($ L ) + 2l R g R cos(A i? - 5 R ) cos($ R ) + g\ + g\) 

where ^ is the CP-eigenvalue of V. For V = {p, co, 0, K*(Ks7i )} the eigenvalue is ^ = 1 
and for V = K*(K L ir°) it is £ = -1. 
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